A MATLAB method of lines template for transport equations
نویسندگان
چکیده
Many environmental problems involve diffusion and convection processes, which can be described by partial differential equations (PDEs). This paper will describe the development of a MATLAB template that generates a numerical solution to PDEs using the method of lines. The template will be applied to various unsaturated flow problems within soil physics to demonstrate the versatility of the method. In particular, the template will generate solutions for three cases (1) one-dimensional Richards’ equation for vertical infiltration; (2) coupled one-dimensional Richards’ equation and solute transport equation for horizontal water and contaminant flow; and (3) two-dimensional Richard’s equation for unsaturated flow over a complex geometry. Where possible, the results from the template will be compared against analytical solutions to determine the accuracy of the numerical solution. In addition, the paper will provide a discussion on possible extensions to the template and future directions. 2003 Elsevier Ltd. All rights reserved.
منابع مشابه
Corrigendum to "A MATLAB method of lines template for transport equations" [Environ. Model. Softw. 19(2004) 603-614]
متن کامل
Solving reaction–diffusion equations 10 times faster
The most popular numerical method for solving systems of reaction–diffusion equations continues to be a low order finite–difference scheme coupled with low order Euler time stepping. This paper extends previous 1D work and reports experiments that show that with high–order methods one can speed up such simulations for 2D and 3D problems by factors of 10–100. A short Matlab code (2/3D) that can ...
متن کاملSequential Implicit Numerical Scheme for Pollutant and Heat Transport in a Plane-Poiseuille Flow
A sequential implicit numerical scheme is proposed for a system of partial differential equations defining the transport of heat and mass in the channel flow of a variable-viscosity fluid. By adopting the backward difference scheme for time derivative and the central difference scheme for the spatial derivatives, an implicit finite difference scheme is formulated. The variable-coefficient diffu...
متن کاملتخمین عملکرد خطوط و ریسک عایقی ناشی از ضربات مستقیم صاعقه با استفاده از روش ترکیبی جدید مبتنی بر روش مونت کارلو
This paper presents a new method, based on analytical relations and Monte Carlo simulation for calculating the lightning performance of overhead lines caused by direct strokes. The aim of the work is to reduce the computing time together with retaining enough accuracy to estimate the line lightning performance and performing the risk analysis of power networks. In the proposed method, some modi...
متن کاملDevelopment of a Computer Model for the Prediction of Littoral Sediment Transport Rate
In this paper, longshore sediment transport in littoral zones is investigated. For investigation of sediment transport in the nearshore zone, the effects of waves, currents and topographical conditions of coast are considered. Linear wave theory was used for the investigation of the wave behaviour. Governing equations of littoral current are continuity and momentum ones. For calculating concent...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Environmental Modelling and Software
دوره 19 شماره
صفحات -
تاریخ انتشار 2004